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We theoretically consider the temperature and density dependent transport properties of 
semiconductor-based 2D carrier systems within the RPA-Boltzmann transport theory, taking into 
account realistic screened charged impurity scattering in the semiconductor. We derive a leading 
behavior in the transport property, which is exact in the strict 2D approximation and provides a 
zeroth order explanation for the strength of metallicity in various 2D carrier systems. By carefully 
comparing the calculated full nonlinear temperature dependence of electronic resistivity at low tem- 
peratures with the corresponding asymptotic analytic form obtained in the T/Tf — > limit, both 
within the RPA screened charged impurity scattering theory, we critically discuss the applicability 
of the linear temperature dependent correction to the low temperature resistivity in 2D semicon- 
ductor structures. We find quite generally that for charged ionized impurity scattering screened 
by the electronic dielectric function (within RPA or its suitable generalizations including local field 
corrections), the resistivity obeys the asymptotic linear form only in the extreme low temperature 
limit of T/Tf < 0.05. We point out the experimental implications of our findings and discuss in the 
context of the screening theory the relative strengths of metallicity in different 2D systems. 
PACS Number : 71.30.-Hh; 73.40. Qv 



I. INTRODUCTION 

A great deal of attention has recently been focused 
on the temperature dependence of carrier (both elec- 
trons and holes, depending on whether the 2D system 
is n-doped or p-doped - in this paper the terminology 
"electron" or "electronic" generically refers to electrons 
or holes depending on the 2D system being considered) 
resistivity, p(T), at low temperatures (and densities), fol- 
lowing the pioneering experimental report by Kravchcnko 
and collaborators [1] that the measured low temperature 
p{T) shows very strong "metallic" temperature depen- 
dence (i.e. p(T) increasing with T) at some intermedi- 
ate densities (the so-called "metallic" or the 2D "metal" 
phase) eventually making a transition to a strongly in- 
sulating state at low carrier densities (n). (At high elec- 
tron densities p(T) shows weak temperature dependence 
[2] similar to 3D metals.) In particular, p{T) could in- 
crease by as much as a factor of three in 2D Si MOS- 
FETs (for n ~ lO^^cm"^) as temperature changes from 
50 mK to a few K [1,3]. While this metallicity (in this 
paper "metallicity" or ""metallic behavior" will exclu- 
sively signify the unusually strong T-dependence of p{T) 
in the metallic phase above the critical density at which 
the system makes a transition to the manifestly insulat- 
ing phase) is by far the strongest in n-Si MOSFET 2D 
structures, the phenomenon has by now been observed 
(with large quantitative variations in the strength of the 
metallicity) in essentially all the existing low density 2D 
semiconductor systems [4-11] such as p-GaAs, p- and n- 
SiGe, Si on sapphire (SOS), n-GaAs, n-AlAs, etc. Our 
work (we concentrate here on n-Si MOS, p-GaAs, and 
n-GaAs, as representative 2D systems) presented in this 
paper deals with the currently controversial issue of un- 



derstanding this metallicity from a theoretical perspec- 
tive. In particular, we use a conventional Fermi liquid 
theory approach in explaining the strong temperature de- 
pendence of p{T) in the metallic phase. We use the well- 
established RPA-Boltzmann transport theory for calcu- 
lating p{T) for 2D carrier systems, taking into account re- 
sistive scattering of the carriers by RPA-screened charged 
impurity random potential. The basic physical picture is 
that of a strongly temperature dependent effective disor- 
der seen by the 2D electrons at low carrier densities due 
to the temperature dependent screening of charged im- 
purity scattering which gives rise to the dominant resis- 
tive mechanism in semiconductors at low temperatures. 
We have earlier obtained qualitative agreement with ex- 
perimental low density p{T) measured in Si-MOS [12], 
p-GaAs [13], SiGe [14], Si-MOS with substrate bias [15], 
and n-GaAs 2D structures [16] using this microscopic 
screening theory approach. For higher carrier densities, 
however, this screening theory is known [18] to provide 
an excellent quantitative description of 2D carrier trans- 
port. 

In this paper we consider, motivated by recent theo- 
retical and experimental development, the leading-order 
temperature dependence of the 2D "metallic" resistiv- 
ity in the low temperature, T/Tp —f 0, limit (where 
Tp — Ep/ks oc n is the 2D Fermi temperature), and pro- 
vide a qualitative explanation for the relative strength of 
metallicity in various 2D systems. Such a unifying qual- 
itative explanation for the relative metallicity strengths 
in different materials has so far been lacking in the liter- 
ature. The 2D "metallic" phase is unusual in the sense 
that the usual 3D metals do not exhibit strongly tem- 
perature dependent resistivity (unless of course there is 
a superconducting transition, not relevant for our con- 
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sideration) at low temperatures (T < 5K), the so-called 
Bloch-Griineisen regime where acoustic phonon scatter- 
ing (the main mechanism contributing to the tempera- 
ture dependence of resistivity in bulk metals) essentially 
freezes out. Although phonon scattering plays a subtle 
(albeit secondary) role [13] in the metallic behavior of 
GaAs-based (both electron and hole) 2D systems, the- 
oretical calculations [17,13,16] definitively show phonon 
scattering to be of little significance in the observed low 
temperature (< IK) 2D metallicity - in fact, in Si MOS- 
based 2D electron systems, where the 2D metallicity is 
most pronounced (and first observed), phonon scattering 
plays no roles whatsoever in p{T) for the experimentally 
relevant regime of T < 5K. Phonon scattering effects, 
which we have considered elsewhere [17,13,16] in provid- 
ing an explanation for the observed non-monotonicity in 
p{T) at intermediate temperatures (T ~ 1 — 5K) in 2D 
n- and p-GaAs systems, are not included in the current 
work since the focus of this paper is the behavior of p{T) 
as T/Tp where phonons surely play no roles. We 
consider only disorder scattering due to random charged 
impurities (and surface roughness scattering, cf. see Sec. 
V) in this work. 

The strong temperature dependent "metallic" resistiv- 
ity in low density 2D systems arise, in our view, from an 
interplay in the disorder scattering between finite tem- 
perature (or, even "high" temperature in the sense that 
T/Tp ~ is not necessarily small as it is in 3D metals 
and could actually be of order unity in low density 2D 
systems for T IK) and density dependent 2D screen- 
ing properties as reflected in the dimensionless parameter 
qTr/'^kp ~ nT^/"^ where qtf and kp are respectively the 
2D Thomas-Fermi screening wave vector and the Fermi 
wave vector [18]. The fundamental difference between 
the semiconductor-based 2D "metals" of interest to us 
in this paper and the usual 3D metals is the great dis- 
crepancy in the magnitudes oiT/Tp and qpF/'i-kp in the 
two systems : In 3D metals T/Tp ~ 10"'' for T ^ IK 
whereas T/Tp ~ 0.1 — 1 in 2D semiconductor systems, 
and qTp/2kp 1 in 3D metals whereas qxp/2kp varies 
between 0.1 and 20 as carrier density is changed in the 
2D systems. In addition static screening has qualita- 
tively different wave vector dependence in 2D and 3D 
systems, leading to the observed strong metallicity in var- 
ious 2D systems. In 3D metals the low temperature re- 
sistivity arises almost entirely from temperature indepen- 
dent short range disorder scattering, which leads to ex- 
ponentially suppressed, 0(e~^/^*'), temperature depen- 
dence in the resistivity, and any residual small tempera- 
ture dependence in p{T) is contributed by phonon scat- 
tering (which produces the well-known Bloch-Griineisen 
behavior, p{T) w po + AT^, where the temperature- 
independent contribution pq arises from short-range dis- 
order scattering whereas the very weak temperature de- 
pendence characterized by the second term arises from 
highly suppressed phonon scattering at low tempera- 



tures). By contrast, low temperature transport (neglect- 
ing weak localization effects) in 2D metallic systems of 
interest to us is dominated mostly by screened disordered 
scattering [i.e., piT) = po + Ap{T) with both po and Ap 
being determined essentially by disorder for T < 5K], 
which can be strongly temperature dependent at low 
densities by virtue of large possible values of the rele- 
vant parameters T/Tp (~ 1) and qTp/2kp (~ 10 - 20) 
at low densities and temperatures in 2D semiconductor 
structures. All localization (as well as interaction effects 
beyond RPA) effects are ignored [2] in this paper. 

The paper is organized as follows. In section II 
we discuss a scaling property (see Appendix A for the 
theory) of the Boltzmann theory resistivity within the 
RPA screened charged impurity scattering model, which 
provides a zeroth order qualitative explanation for the 
strength of the transport metallicity in 2D carrier sys- 
tems. In section III we discuss the asymptotic low tem- 
perature {T/Tp 0) behavior of the calculated resis- 
tivity in the RPA-Boltzmann model comparing it quan- 
titatively with the full temperature dependent resistiv- 
ity (in the same modal) in order to estimate the regime 
of validity of the leading-order temperature expansion 
of resistivity. In IV we consider the various solid state 
physics effects (e.g. the quasi-2T) nature of the semi- 
conductor layer, the long-range or the short-range na- 
ture of the bare scattering potential) on the 2D transport 
properties. In section V we provide a critical compari- 
son between our realistic (but theoretically approximate) 
RPA-Boltzmann 2D transport theory and a set of re- 
cent experimental results in Si inversion layer, concluding 
that our theory, without any adjustment of parameters 
and/or ad hoc theoretical refinement, describes well the 
observed experimental temperature dependence down to 
a carrier density of about 5 x IG^^cm"^ — for lower den- 
sities the agreement between experiment and theory is at 
best qualitative with the actual temperature dependence 
of p{T) being stronger than the calculated p{T). We con- 
clude in VI with a discussion of the implications of our 
results. 

II. DENSITY-TEMPERATURE (go, t) SCALING 
OF METALLICITY 

In the RPA-Boltzmann theory (cf. Appendix A) the 
dimensionless Vs parameter (the so-called Wigner-Seitz 
radius) characterizing the electron-electron interaction 
strength in the 2D system does not play a fundamen- 
tal role in determining the temperature dependence of 
p{T) except so far as determines the dimensionless 
parameters T/Tp (= t) and qrp/^kp (= go) through the 
carrier density. We believe that the fundamental minimal 
parameters determining the zeroth order 2D metallicity, 
i.e. the temperature dependence of p{T,n) in the pu- 
tative metallic phase, are t and qo- In particular, for 
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2D systems it is easy to show that qo = Vg/V^ and 
t = {ksT/ Ry){gv/'i)r1, where gv is the vahey degen- 
eracy of the relevant semiconductor material {gv = 2 
for Si(100)-MOS structures, g^ = 1 for p- and n-GaAs 
systems), = {im)^^/'^ / as is the usual dimension- 
less Wigner-Seitz density (or interaction) parameter with 
qb = nf? jrae? as the effective semiconductor Bohr ra- 
dius (k and m are the background dielectric constant 
and the carrier effective mass respectively), and Ry = 
I (2KaB) is the natural atomic energy unit (effective Ry- 
dberg) for the semiconductor. We find that the existing 
experimental data for the metallicity in various 2D semi- 
conductor systems approximately obeys the 2-parameter 
scaling behavior Ap{T,n)/ po ~ F{t,qo), where F is a 
smooth and approximately universal function of go and t 
for all 2D "metallic" systems, as implied by the screening 
theory. In particular, a direct consequence of this theo- 
retical prediction is that the temperature dependence of 
p(T) should correlate with the parameter go in different 
materials when expressed in terms of the dimensionless 
temperature t. (We note that the functional dependence 
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of qa ^ g'v is not only different from the dimension- 
less density parameter r^, but also from the dimensionless 
ratio of the Fermi energy to the Coulomb energy in the 
system which goes as gv^s-) Theoretical details and the 
equations for our 2D RPA-Boltzmann theory are given in 
the Appendix A of this paper where we derive this scaling 
law. We note that the scaling behavior Ap/po ^ F{t, qo), 
with t = T/Tp and go = QTp/'^'kF, derived in this paper 
is exact for RPA screened charged impurity scattering 
in the ideal 2D limit in contrast to various other scal- 
ing behaviors (e.g. refs. [12,13,15]) discussed earlier in 
the literature which are all approximate scaling behavior 
valid only in limited range of parameters. 

This zeroth order functional dependence on t and go, 
Ap{T,n) = Ap{t,qo), in fact, provides a minimal expla- 
nation for the observed strong variation in the 2D metal- 
licity not only for various densities in the same material 
but for different materials at equivalent densities - for 
example, Si-MOSFET based 2D electron systems man- 
ifest much stronger metallicity compared with ^GaAs 
or n-GaAs based 2D systems even at the same value 
because g„ = 2 (1) in Si (GaAs). Taking to be the 
critical parameter determining metallicity fails to explain 
why the 2D p-GaAs system shows weaker metallic behav- 
ior, even though it typically has much larger values 
(fg ~ 15 — 40) than the Si system. The screening the- 
ory provides a simple explanation for this observation by 
virtue of go being larger (by a largo factor of VS) in Si 
than in GaAs for the same value. Thus, Si MOSFET 
2D metallicity at = 10 should be approximately com- 
parable to a GaAs 2D metallicity at rg = lOVS ~ 28 
when expressed as a function of the dimensionless tem- 
perature variable t = T/Tp. This predicted correspon- 
dence in the relative metallicity in terms of go and t is 



consistent with the experimental observations in these 
systems. We emphasize, however, that this correspon- 
dence is expected to work only on a zeroth-order quali- 
tative level and should not, by any means, be construed 
as a precise quantitative prediction. This is discussed 
with theoretical details in the appendix. For example, 
the form factor effects associated with the quasi-2D sub- 
band quantization do not scale with the density parame- 
ter go, and will necessarily affect different systems in dif- 
ferent manners since the associated effective masses, the 
dielectric constants, the depletion charge densities, and 
the confinement potentials, which together determine the 
form-factor, are different in different systems. Similarly, 
the bare (i.e., unscreened) disorder will certainly de- 
pend on the system varying qualitatively among different 
systems and materials, which could lead to substantial 
quantitative deviations from our predicted p(go, t) scaling 
based on the temperature dependent screening argument. 
For example, in Si MOSFETs transport is dominated by 
interface scattering - both by long-range potential scat- 
tering due to ionized impurities located near the Si-Si02 
interface and short-range scattering by interface rough- 
ness fluctuations inherent at the Si-Si02 interface. In 
GaAs structures (both p and n), scattering by (uninten- 
tional) background charge impurities (and less so by re- 
mote charge dopants) dominates at the highest mobilities 
with phonon scattering being non-negligible (although 
small) down to T < IK . In addition, p-GaAs struc- 
tures most likely also have significant spin inter-subband 
scattering within the spin-split valence bands. The non- 
universal quantitative effects associated with subband 
quantization and disorder potential make the simple 2- 
parameter screening picture of Ap{T, n) dependent only 
on go and t quantitatively inaccurate, but the simple pic- 
ture applies surprisingly well on a zeroth-oder qualita- 
tive level as can be verified by comparing the experimen- 
tally observed metallicity strengths in Si MOS, p-GaAs, 
and n-GaAs structures where the metallicity scales ap- 
proximately with go and t provided T is low enough so 
that phonon effects could be ignored in the GaAs sys- 
tem. The crucial new point we are making in this paper 
is that strong metallicity manifests itself in low density 
2D systems because the control parameters go oc n~^/^ 
and t oc n^^ are large only for low carrier densities and 
not because rg is large - for example, n-Si MOS system 
and n-GaAs system show more than an order of magni- 
tude different metallicities for the same value because 
go and t are much larger in Si than in n-GaAs {Ap/ po in- 
creases by a factor of 3 in Si MOS systems of refs. [1] and 
[15] for r<, > 10 whereas it increases only by about 25% 
in the n-GaAs system of refs. [16] for > 10). Within 
a specific materials system, however, the go dependence 
of the resistivity becomes completely equivalent to an rg 
dependence (since gi, is a constant for a given system, 
and go = gt^^rs/V^) as one would expect — it is only in 
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FIG. 1. Calculated Ap — p(T,n) — po, where 
po = p{T = 0), for various qo = qTv/'^kp — 5, 10, 15, 20 (from 
the bottom) as a function of t = T/Tf for Si MOSFET. The 
solid (dashed) lines indicate resistivities from interface (bulk) 
charged impurities. 

comparing different systems (e.g. Si and GaAs) that qo 
and Vs dependence are not equivalent. 

In Fig. 1 we show our calculated p{T,n), depicted 
in terms of the dimensionless variables qo and t, for a 
strictly 2D system using Si MOS parameters (e.g., effec- 
tive mass, dielectric constant, etc.) and assuming the 
disorder scattering to be entirely due to finite tempera- 
ture RPA-screcned charged impurity scattering. Results 
shown in Fig. 1 demonstrate the importance of the di- 
mensionless screening parameter qo (= qTv/^kp) as the 
relevant control parameter in determining the strength 
of metallicity in 2D systems. In particular the maximum 
relative change in the resistivity, Ap/po, ranges from 
about 40% for go = 5 (which corresponds to a Si (100) 
inversion layer carrier density of n = 11 x 10^^ cm~^, 
Tf = SOK) to almost 300 % for qo = 20 (corresponding 
to n = 7 X IQi" cm-^ Tp = 5K) as t = T/Tp changes 
from to 1 in Fig. 1. Therefore, the RPA results of Fig. 
1 indicate an increase of resistivity by about 10% and 
300% respectively for carrier densities 1.1 x 10^'^cm~'^ and 
7x 10i°cm-2 in Si MOSFETs for a change in T ofO-5K, 
assuming that the system remains metallic. It should be 
noted that the maximum in Ap/ po shifts to higher val- 
ues of t (= T/Tp) for higher (lower) values of qTp/2kp 
(n), and phonon effects (ignored in our consideration) 
will play increasingly important quantitative role in the 
temperature dependent resistivity for T > 5K. These 
two facts together make the metallic behavior relatively 
even more important at lower densities, or equivalently, 
higher values of qo- 

A comparison between solid and dashed lines in Fig. 1 
shows the qiiantitative importance of the nature of im- 
purity scattering in determining the temperature depen- 
dence of resistivity: in general, charged impurities ran- 



domly distributed in the 2D layer itself (bulk disorder 
shown by dashed lines in Fig. 1) lead to stronger tem- 
perature dependence than interface disorder (solid lines) 
associated with charged impurities distributed randomly 
at the semiconductor-insulator interface. This is pre- 
cisely what is expected since screening effects should be 
the strongest when charged impurities and the carriers 
are distributed in the same region of space with no spa- 
tial separation. It may also be worthwhile to mention in 
the context of Fig. 1 that the experimental Si inversion 
layer systems [1,3] manifesting the most dramatic metal- 
licity {i.e. large changes in Ap/po as a function of tem- 
perature) typically have n < 10^^cm~'^ corresponding to 
qo 15 — 20 in Fig. 1, thereby showing a relative tem- 
perature dependent change in resistivity of about 100 — 
300 % as T/Tp varies from zero to 0.5. Thus the results 
in Fig. 1 are in reasonable qualitative agreement with ex- 
perimental results as we have emphasized elsewhere [12]. 
Scattering mechanism not included in the theory (e.g. 
surface roughness scattering) and higher-order interac- 
tion effects will certainly modify the quantitative details 
of the results, but it is gratifying to see that a zeroth- 
order Boltzmann transport theory including only RPA 
screened charged impurity scattering provides a reason- 
able qualitative description of the observed metallicity. 

In Fig. 2 we show the calculated Ap/ po for three dif- 
ferent 2D systems for a comparison of their metallicity: 
(100) n-Si inversion layer, p-GaAs heterostructure, and 
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FIG. 2. The main figure shows the scaled resistivity as a 
function of t = T/Tp for Si-MOSFET, p-GaAs, and n-GaAs 
with a fixed qo = 15. For pure 2D system the calculated 
resistivities yield perfect scaling for all systems (thin solid 
line). Thick lines show the scaling for quasi-2D systems. In 
inset (a) the comparison of the metallicity for fixed density 
n = 10^ ^cm~^ is given as a function of a temperature. In 
inset (b) we show the change of resistivity for fixed value of 

= 8 as a function of t = T/Tp. 
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n-GaAs heterostructurc. Wc have shown the results for 
the same values of carrier density, parameter, and go 
(= iTF/^kp) parameter in order to emphasize the impor- 
tance of the dimcnsionlcss screening parameter qo in de- 
termining the temperature dependence. As emphasized 
above, for the same values of carrier density (see inset 
(a)), the mctallicity is the strongest (weakest) in the Si 
(n-GaAs) system with the p-GaAs system being interme- 
diate. For the same value of Vg parameter (see inset (b)) 
the metallicity depends strongly on the valley degener- 
acy, but weakly on the other parameters. The mass of 
the carrier dose not change the metallicity for the same 
parameter, which gives exactly the same resistivity be- 
havior for p-GaAs and n-GaAs if we use the same mate- 
rial parameters except mass. In the main figure we show 
the metallicity for equivalent go = 15 parameter values 
in all three systems. The thin solid line shows that the 
temperature dependence is scaled perfectly for the pure 
2D systems, as proven in the appendix. For quasi-2D 
systems, however, the scaled resistivity expressed as a 
function of dimensionless temperature t = T/Tp is ap- 
proximately similar in the three systems without per- 
fect scaling. The deviation of scaling is mostly due to 
form factor effects. Since the electron effective mass in 
n-GaAs is small (m = 0.067rne), extremely low values 
of carrier density are required in 2D n-GaAs systems for 
observing appreciable screening-induced temperature de- 
pendence as has been recently reported in Ref. [16]. 

III. LOW TEMPERATURE ASYMPTOTIC 
BEHAVIOR 

Since the behavior of Ap{T,n) for finite T is neces- 
sarily non-universal for reasons discussed above, recent 
attention has focused on the very low temperature be- 
havior of p{T) in the t = T/Tp limit. In particular, 
it was realized a long time ago [19,20] that Ap{T, n) de- 
rived from the Boltzmann theory within the relaxation 
time approximation for RPA screened disorder scattering 
has the following expansion in 2D systems: 

Ap(r, n)/po ~ Cit + Cz/2t"^ + ■■; (1) 

in the asymptotic t ^ Q regime. The coefficients Ci and 
C3/2 were first (but incorrectly) calculated by Gold and 
Dolgopolov [20]. We have re-calculated C\ and C3/2, and 
find the original results in Ref. [20] to be incorrect. We 
calculate the correct coefficients (for RPA screening) to 
be: 

Ci=2{l + llqofr\ (2a) 
6-3/2 = 2.646 (1 + 1/go/)"', (2b) 

where go = qTE/'^kp (as defined above) and / = f{2kp) 
is the appropriate quasi-2D subband form-factor at the 



wave vector 2kF (in general / < 1 with the strictly 2D 
limit being / = 1, see appendix A). Our calculated Ci 
agrees (in the strictly 2D limit of / = 1) with the re- 
cent (Hartree) result given in Ref. [21] and disagree with 
that of Ref. [20] whose incorrectly calculated C\ is larger 
by a factor of 2 In 2 (i.e. about 40% higher). Our cal- 
culated C3/2 ~ 2.65 is 28% smaller in magnitude than 
the incorrect value {~ 3.4) quoted in Ref. [20]. We note 
that the errors in Ref. [20] led to a large overestimate of 
the asymptotic temperature dependence of p{T) in the 
Gold-Dolgopolov formula [20], which we now correct. An 
important recent theoretical development in the subject 
has been the demonstration by Zala et al. [21] that the 
leading-order linear result given in Eq. 1, in fact, survives 
(albeit with C\ replaced by an unknown Fermi liquid pa- 
rameter) inclusion of higher-order electron-electron inter- 
action terms in the theory, of which screening is only one 
particular aspect. In particular, the leading-order tem- 
perature dependence in the interaction theory of Zala et 
al. contains the C\ term of our Eqs. (1) and (2) as the 
so-called Hartree term in the language of Ref. [21]. 

A thorough understanding of this first order linear 
term in the theory has taken on significance in view 
of the existence of the Zala et al. [21] work, and even 
more importantly, because of the several recent attempts 
[11,22,23] to compare experimental results to the interac- 
tion theory. The interaction theory considerably extends 
the screening theory (through the inclusion of higher- 
order interaction corrections), but is unfortunately con- 
strained at this stage to only the leading order result 
in t, and therefore applies only at very high (low) den- 
sities (temperature) so that the constraints f <C 1 and 
Ap/po <C 1 are satisfied. The two theories are thus com- 
plementary - the screening theory being an approximate 
theory (because it includes only the screening part of the 
electron-electron interaction) for all t (in fact, its accu- 
racy improves with increasing t since the RPA becomes 
exact in the classical high-temperature limit) whereas the 
interaction theory is presumably an exact leading-order 
in t (within the perturbative Landau Fermi liquid theory 
scheme) theory as i ^ 0. This obvious complementar- 
ity of these two approaches has not been emphasized in 
the recent literature where some recent publications have 
even presented the misleading and incorrect viewpoint of 
these two approaches as mutually exclusive competing 
theories. It is important to emphasize that the inter- 
action theory [21], by construction, applies only when 
the temperature correction to the T = conductivity is 
small, i.e. the theory of ref. [21] is a leading order theory 
for small temperature corrections to the T = conduc- 
tivity as T/Tp — > 0. By definition, therefore, this in- 
teraction theory cannot explain the strong metallicity or 
the large temperature-dependent changes in the resistiv- 
ity reported in the literature. The real significance of ref. 
[21] is strictly theoretical — it shows that the "metallic" 
behavior given by Eq. (1) within the RPA-Boltzmann 
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theory survives higher-order electron-electron interaction 
corrections in the T/Tp limit (provided that the 
"ballistic" transport condition ksT h/r, where r is 
the T = transport relaxation time (t = m/ne^po), 
is satisfied, i.e., the temperature is in the intermediate 
range Tp » T /i/rfes). We note that the necessary 
condition minimally required for a comparison between 
experiment and the interaction theory is that (1) the 
experimental conductivity must show a linear temper- 
ature correction in the intermediate temperature regime 
Tp ^ T ^ h/rkB, and (2) the temperature correction to 
the T = conductivity must be small. Most experiments 
on 2D transport do not satisfy these necessary conditions 
for comparison with the interaction theory, most particu- 
larly because the measured conductivity essentially never 
manifests a linear temperature dependence except at very 
high densities where the RPA-Boltzmann semiclassical 
transport theory is quantitatively accurate. 

In view of the complementary nature of these two the- 
ories it becomes important to ask about the regime of 
validity of the linear approximation inherent in the in- 
teraction theory. This issue becomes particularly crucial 
since most of the existing p{T) data in the putative 2D 
metallic phase does not follow a linear temperature de- 
pendence over any appreciable temperature regime in the 
lowest temperature range (i.e. T/Tp <ti 1). The situation 
becomes more complicated with the realization that the 
'strong' condition (the 'weak' condition being T <^Tp) 
for the validity of the interaction (as well as the screen- 
ing) theory is that Td <^T <^Tp, where Tjj ^ h/{TkB) 
is roughly the so-called Dingle temperature with r be- 
ing the T = limit of the transport relaxation time [i.e. 
T = m/{ne'^po)]. In the screening theory Td cuts off 
the temperature dependence of screening for T < To 
(making the disorder to be effectively temperature inde- 
pendent for T <^ Td) leading to a suppression in the 
temperature dependence of p{T) whereas the interaction 
theory is by construction a theory of "ballistic" trans- 
port developed in the h/irkBT) ^ 1 limit (and then the 
T /Tp <C 1 limit is taken to obtain explicit asymptotic 
results). The T^-cutoff in the temperature dependence 
of p{T) for T < is extremely well-motivated physi- 
cally within the screening theory, and has been discussed 
in details in the literature [12,13,19] as the reason for the 
need of low disorder (or equivalently, high mobility with 
concomitant low values of Tpt) samples to observe 2D 
metallicity (low density is also required in the screening 
theory for strong metallicity so as to produce large values 
of T/Tp at low temperatures and to make qxp/^kp large 
enough to have strong screening effect). We will mostly 
ignore the Dingle temperature effects in the theory by as- 
suming Td « 0, but in comparing experimental data to 
the interaction theory the ballistic limit is an important 
constrain to remember. In particular, in many exper- 
imental situations the constraint /i/(rfcB) <C T <C Tp 
necessary for the application of the ballistic limit inter- 




FIG. 3. p(T) for a Si-MOSFET with density n =56.0, 12.0, 
7.0 X 10^° cm~^ (from the bottom). Solid lines indicate the full 
RPA screening theory, dashed lines the asymptotic approxi- 
mations with only the linear term, and dot-dashed lines with 

both the 0(r) and 0{T^'^) terms. For clarity wo use offset by 
0.5 for lower densities. The inset shows p{T) on an expanded 
temperature scale. 

action theory may not even exist [23]. 

To address the important question of the regime of 
validity of the asymptotic linear approximation we have 
carried out a careful comparison between the linear-T 
approximation (Eqs. 1 and 2) and the full numerical 
calculation of p{T) within the RPA-Boltzmann theory. 
(These results are presented in Figs. 3 and 4.) The im- 
portant conclusion drawn from this comparison is that 
the regime of validity of the asymptotic linear formula, 
Ap/po ~ 0{T/Tf), is extremely restricted, and at least 
for the RPA screened disorder scattering, the linear ap- 
proximation holds only in the very low temperature limit 
of T/Tp < 0.05, which, in general, is not in the balhstic 
regime except at very high densities where the metallic- 
ity is very weak. We therefore conclude that the low- 
temperature asymptotic linear regime, while being of 
substantial theoretical interest, is not of much experimen- 
tal relevance since Ap ^ T only for temperatures (densi- 
ties) much too low (high) to be of experimental interest. 
Our conclusion is based on a comparison of the asymp- 
totic analytic linear-T formula with the full T-dependent 
calculation only for the RPA screening theory (because 
this is the only theory where both the asymptotic result 
and the full result can be calculated), but we believe our 
conclusion to be quite generally valid, and even for the 
interaction theory we expect the low-T regime of validity 
of the linear T formula to be too restricted to be of much 
experimental relevance. This is consistent with the exist- 
ing experimental results where the conductivity is never 
precisely linear (at low temperatures) essentially at any 
carrier densities, except perhaps at very high carrier den- 
sities where the RPA-Boltzmann theory should 
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FIG. 4. Calculated p{T) as a function of t = T/Tf (a) 
for p-GaAs with density n =21.8, 5.45, 1. 36x10^° cm~^ (from 
the bottom) and (b) for n-GaAs with density n =6.1, 1.5, 
0. 38x10^ cm~^ (from the bottom) by considering only scat- 
tering from a uniform random distribution of ionized impu- 
rities in the quantum well. Solid lines indicate the full RPA 
screening results, dashed lines the asymptotic approximations 
with only the linear term, and dot-dashed lines with both the 
0{T) and 0(T^''^) terms. Offset by 0.5 for lower densities is 
used. The inset in (a) shows p{T) for charged impurities at 
interface. The inset in (b) shows p{T) for scattering by re- 
mote charged dopants at d = lOOA (solid lines) and d = 300A 
(dashed lines) on an expanded temperature scale. 

be quantitatively valid. 

In Fig. 3 we show our calculated p{T) for three den- 
sities in the (lOO)Si MOS 2D electron system in both 
the full RPA screening theory and in the asymptotic ap- 



proximations keeping only the leading-order linear term 
and both the 0{T) and ©(T^/^) terms in Eq. (1). We 
have assumed scattering from a random distribution of 
charged impurities located at the interface, and subband 
form-factor effects have been included in both calcula- 
tions equivalently. The inset shows p{T) on an expanded 
temperature scale. 

In Fig. 4 we show similar comparisons between the 
calculated full p{T) and the asymptotic analytic approx- 
imations for 2D GaAs holes and electrons in their exper- 
imentally relevant density regimes of interest. For these 
two systems we show results for two different impurity 
scattering mechanisms for the sake of completeness (also 
for the sake of consistency with the experimental GaAs 
systems where interface charged impurity scattering is 
typically less important than in Si MOSFETs). In par- 
ticular, we give results for scattering by a uniform ran- 
dom distribution of (unintentional) background ionized 
impurities (which are invariably present and are usually 
the dominant scattering centers in the GaAs samples of 
the highest mobilities) and by remote charged dopants 
(assumed to be randomly distributed in a 2D plane sep- 
arated a modulation distance of d from the 2D electron 
layer), the so-called modulation doped situation. It is ob- 
vious from Fig. 4 that 2D n-GaAs system not only has 
the weakest temperature dependence but also exhibits 
essentially no clear linear temperature regime. This is 
consistent with the very weak screening properties (and 
large values of Ep) of 2D n-GaAs because of its very 
low band effective mass. The real importance of the ac- 
tual random impurity distribution in the system (usually 
not known and has to be inferred from a comparison of 
the transport data with theoretical calculations assuming 
specific impurity distributions) in affecting (both quali- 
tatively and quantitatively) the p(T) behavior is appar- 
ent in the results of Fig. 4: The strongest temperature 
dependence arises in the situation where the charged im- 
purities are randomly distributed in the 2D layer of the 
carriers and the weakest T-dependence arises in the mod- 
ulation doped situation (particularly for densities low 
enough so that 2k pd > 1) with the interface disorder 
case being intermediate. This dependence on the details 
of disorder is easily understood within the screening the- 
ory by considering the role of 2kp scattering in transport; 
Modulation doping with 2kpd > 1 essentially completely 
suppresses large momentum 2kp scattering even at low 
temperatures because of the e^^/cFrf term in the form fac- 
tor drastically reducing the temperature dependence due 
to screened impurity scattering whereas charged impuri- 
ties distributed randomly in the 2D layer maximizes 2kp 
scattering for a given system. The interaction theory 
of Zala et al. [21] had to make the drastic approxima- 
tion of a zero-range white-noise impurity disorder po- 
tential, thus drastically (and artificially) enhancing the 
2fcir-scattering. In real systems, other things being equal 
(i.e. the mobility, the density, and the 2D system), there 
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T/Tp 

FIG. 5. The calculated pure-2D (dashed lines) and 
quasi-2D (solid lines) temperature-dependent resistivity for 
screened bulk charge impurity scattering (a) for Si-MOSFET 
with n = 10, 20x 10^° cm"^, (b) for p-GaAs with p = 
2.0, S.OxlO^^cm-^ and (c) for n-GaAs with n = 0.5, 
LOxlO^^cm"^. The lower lines represent higher densities. 

would be a strong dependence of the detailed behavior of 
p{T) on the actual random impurity distribution in the 
system through the form factor effect. (Note that the 
low-T asymptotic linear formula does not depend on the 
impurity distribution or on the range of disorder scatter- 
ing, but its temperature regime of validity may very well 
depend on the nature of scattering in the 2D system.) 

IV. SOLID STATE PHYSICS EFFECTS: 
QUASI-2D LAYER WIDTH AND BARE 
IMPURITY POTENTIAL RANGE 

We now discuss the quantitative significance of various 
solid state physics effects on the 2D metallicity. The spe- 
cific effects we discuss are the quasi-2D nature of the 
semiconductor layers under consideration and the na- 
ture of bare impurity disorder {i.e. long range versus 
short range) in determining the temperature dependence 
of 2D resistivity. (Both of these effects are ignored in 
Ref. [21]. The inclusion of the quasi-2D form- factor ef- 
fect in the interaction theory is straightforward, but the 
inclusion of long-ranged bare disorder, e.g. charged im- 
purity scattering, in the interaction theory is nontrivially 
difficult.) In Fig. 5 we compare the calculated 2D and 
quasi-2D temperature-dependent resistivity for screened 
bulk charge impurity scattering, finding that the metal- 




FIG. 6. The calculated temperature dependent resitivi- 
ties for realistic long-range charged impurity scattering (solid 
lines) and for hypothetical short-range delta-function impu- 
rity scattering (dashed lines) (a) for Si-MOSFET with n = 10, 
20x 10^° cm^^ (b) for p-GaAs with p = 2.0, 5.0xlO^Vm-^ 
and (c) for n-GaAs with n = 0.5, l.Ox 10^° cm"^. The lower 
lines represent higher densities. 

licity is substantially overestimated in the strictly two- 
dimensional approximation. Thus, it is essential to in- 
clude the quasi-2D nature of the semiconductor system 
in order to obtain quantitatively accurate temperature 
dependence although the 2D approximation is qualita- 
tively correct. 

In Fig. 6 we show a comparison between the results cal- 
culated for realistic long-range charged impurity scatter- 
ing and hypothetical short-range delta-function (in real 
space) impurity scattering (both equivalently screened 
by the finite temperature RPA dielectric function of the 
2D carriers). We point out that the dominant disor- 
der in 2D semicondTictor systems arises from long-range 
(Coulomb) charged impurity scattering. In Fig. 6, the 
calculated temperature dependence is obviously substan- 
tially stronger for short-range bare disorder. It is impor- 
tant to emphasize in this context that the leading-order 
asymptotic dependence of Ap/po on T/Tp is independent 
of the range of the bare disorder since it depends (see 
Eqs. 1 and 2) on the electron-impurity interaction only 
through Vc-i{2ki ), i.e. through the constant (momen- 
tum space) impurity potential defined at the wave vector 
q = 2kF, making long- and short-range bare disorder 
completely equivalent for the low temperature asymp- 
totic temperature dependence. On the other hand, the 
full temperature dependence depicted in our Fig. 6 man- 
ifestly demonstrates that, except at the lowest values of 
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t = T /Tp < 0.1, the actual temperature dependent resis- 
tivity depends quite significantly on the range of bare dis- 
order with short-range bare disorder providing substan- 
tially stronger temperature dependence than the long- 
range Coulomb disorder due to random charged impu- 
rity scattering. Since we have already argued in Sec. Ill 
that the asymptotic linear temperature dependence does 
not really apply, except at extremely low temperatures 
(or high carrier densities), the range of impurity disorder 
potential takes on special significance in the theoretical 
analyses. 



V. COMPARISON WITH EXPERIMENT 

Although the primary goal of this article is to estab- 
lish and clarify certain theoretical principles (e.g. the va- 
lidity of the asymptotic low temperature expansion, the 
strength of metallicity as reflected in the qq, to scaling 
behavior, the importance of various realistic solid state 
physics effects such as the long-range versus the short- 
range nature of the bare impurity disorder potential or 
the quasi-2D nature of the semiconductor systems of ex- 
perimental interest) within the RPA-Boltzmann theory 
of 2D carrier transport in semiconductor systems, it is 
important to ask about the empirical validity of the ze- 
roth order RPA-Boltzmann theory in the context of the 
large amount of the available temperature dependent 2D 
transport experimental data. We have in fact carried out 
a number of comparisons between our theory and the 
experimentally measured temperature dependent resis- 
tivity in several different 2D systems of current interest 
[12-16,23,25-27]. It is important in this context to re- 
member that such comparisons between experiment and 
theory is necessarily qualitative in nature since the ex- 
perimental data, even in the same material system, show 
strong sample to sample variations, and the measured 
p{T, n) in different Si MOS systems (or for that mat- 
ter, in different n-GaAs or p-GaAs systems) often have 
significantly different quantitative (and sometimes even 
qualitative) dependence on temperature and carrier den- 
sity. This is understandable since different samples in 
general may have significantly different bare disorder po- 
tential and system parameters, and therefore a univer- 
sal quantitative behavior of p(T, n) cannot be expected 
(since there is no universal quantitative behavior for the 
experimental data themselves). 

The important question therefore is the extent to 
which a particular theory explains the qualitative behav- 
ior of the observed p{T, n) in 2D systems. The RPA- 
Boltzmann theory discussed in this paper is unique in 
this respect since it is the only theory which is capable 
of producing the full temperature and density dependent 
2D resistivity which can, in principle, be compared with 
the experimental data. The theory is still not uniquely 
defined since all the system parameters (e.g. effective 
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FIG. 7. (a) Calculated p(T, n) within the RPA-Boltzmann 
transport theory assuming only interface charged impurity 
scattering, for all the densities corresponding to Fig. 1(a) 
(Si-22 sample) in ref. 22. (b) Calculated p{T,n) with 
long-range Coulomb impurities randomly distributed at the 
interface (dashed lines) and the ideal zero-thickness 2D layer 
approximation with zero-range 5-function potential bare im- 
purity scatterers (solid lines) for six representative carrier den- 



sities n = 5.7, 6.9, 8.7, 11.7, 
bottom) . 



16.5, 35.7xl0"'cm-^ (top to 



mass, depletion charge density, the precise metallic car- 
rier density, etc.) are never really accurately known, and 
more importantly, the detailed quantitative parameteri- 
zation of the bare disorder potential is never available for 
any sample. In general, there could be many indepen- 
dent sources for the bare disorder (e.g. bulk and inter- 
face charged impurities, surface roughness, remote impu- 
rities, intervalley scattering, phonons, alloy scattering), 
and with a sufficient number of adjustable free param- 
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eters characterizing different kinds of disorder assumed 
to exist in a system one may very well be able to obtain 
essentially complete (but, meaningless) agreement with 
experimental data for any particular sample. 

With all these caveats in mind we carry out a compar- 
ison between our theory and a recent set of experimental 
data [22] for p(T, n) in a Si MOS system. The specific 
sample we choose for this comparison is the Si- 2 2 sample 
(see Fig. 1(a) of ref. [22]) with a quoted "peak mobility" 
of 33,000 cm?' /Vs. In ref. [22] a scries of p{T) curves, for 
T K. O.SJsT — 5.0-ftr, are presented for this sample in the 
2D carrier density range of n = 5.7 — 35.7 x lO^^cm"^ 
with the typical resistivity values spanning between 0.002 
h/e^ to 0.02 h/e^. The temperature induced fractional 
change in p{T) in this 0.5 — S.OJf temperature window 
spans between a few percent at higher densities to about 
35% at the lowest densities. Visually the p(T) curves of 
Fig. 1(a) in ref. [22] corresponding to the Si-22 sample all 
seem to obey the asymptotic linear temperature depen- 
dence although, in general, there is significant deviation 
from linearity both at the lowest (< Q.^K) and at the 
highest (< bK) temperatures. 

In Figs. 7 and 8 we show our calculated p{T,n) for 
the Si-22 sample [22] within the RPA-Boltzmann trans- 
port theory assuming screened disorder scattering. In 
Fig. 7(a) we show the calculated RPA-Boltzmann re- 
sistivity, assuming only interface (or oxide) charged im- 
purity scattering, for all the densities corresponding to 
Fig. 1(a) in ref. [22]. The results shown here look qual- 
itatively similar to the experimental data [cf. Fig 1(a) 
in Ref. [22]] except that the temperature dependence is 
somewhat weaker in the theory. In Fig. 7(b) we show 
the theoretical results, p{T) with T = 0.25 — QK, for six 
representative carrier densities n = 5.7, 6.9, 8.7, 11.7, 
16.5, 35.7xl0^^cm^^ (top to bottom in Fig. 7) using im- 
purity scattering as the only resistive mechanism in the 
system, comparing the quasi 2D realistic situation and 
a long-range bare impurity potential with the idealized 
2D situation and a short-range bare impurity potential. 
The two sets of calculatcid results in Fig. 7 correspond- 
ing to the realistic quasi-2D system (with finite 2D layer 
thickness, semiconductor-insulator dielectric mismatch, 
etc.) with long-range bare Coulomb impurities randomly 
distributed at the interface (dashed lines) and the ideal 
zero-thickness 2D layer approximation with zero-range 
(5-function potential bare impurity scatterers (the model, 
for example, of ref. [21]) providing the resistive scatter- 
ing mechanism (solid lines). We note several important 
features of the results presented in Fig. 7: (1) The theo- 
retical results for p{T) appear to be approximately linear 
in T for all the carrier densities, being qualitatively very 
similar to the experimental data shown in Fig. 1(a) of 
ref. [22]; (2) in spite of this qualitative linearity of p{T), 
the actual curves are nonlinear except at the highest 
densities; (3) an approximate measure (actually a lower 
bound) of this nonlinearity is the difference between the 
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FIG. 8. Calculated p(T) for three carrier densities (a) 
n = 21.3xl0"cm-^ (b) n = 10.5xl0"cm-^ (c) n = 
5.7 X 10^^ cm~^, including realistic surface roughness and in- 
terface (oxide) charged impurity scattering. In (c) the thin 
(thick) solid line represents the result using the band effective 
mass (the measured effective mass from Ref. [24]). Here A is 
the average displacement of the interface and A is of the order 
of the range of its spatial variation in the direction parallel to 
the surface. 



10 



solid and dashed lines in Fig. 7 since the slopes of the 
two lines are the same in the T ^ limit where both 
are strictly linear in T/Tp as discussed above in section 
III; (4) the calculated results, while being qualitatively 
similar to the data [22], disagree quantitatively with ex- 
periment with the quantitative disagreement increasing 
with decreasing carrier density; (5) the unrealistic strictly 
2D approximation with zero-range (5-function bare dis- 
order potential (solid lines in Fig. 7) actually agrees 
much better with the experimental data than the realis- 
tic quasi-2D calculation with charged Coulomb impurity 
scattering as the bare disorder (dashed lines) since the ex- 
perimental data of rcf. [22] shows stronger temperature 
dependence quantitatively than the theoretical results of 
Fig. 7 (again this is consistent with the discussion of Sec. 
IV above where (5-function bare impurity disorder and/or 
strict 2D approximation turns out to give much stronger 
temperature dependence than the realistic quasi-2D sys- 
tem with charged Coiilomb scattering) — we emphasize, 
however, that both the solid and the dashed lines have 
the same slope in the T/Tp — > limit reinforcing the lack 
of linearity in the p{T) behavior. 

To go beyond the qualitative agreement with experi- 
ment depicted in Fig. 7, we consider in Fig. 8 a more 
realistic situation taking into account, in addition to the 
charged Coulombic impurity disorder, the surface rough- 
ness scattering by the Si-Si02 interface which has been 
known for a very long time [18] to be the dominant scat- 
tering mechanism affecting carrier transport in Si MOS- 
FETs at relatively "higher" carrier densities whence the 
2D electron gas resides rather close to the Si-Si02 inter- 
face making the roughness scattering to be significant. 
We include the surface roughness scattering in our the- 
ory in the standard manner [18], screening it within the 
RPA theory similar to the charged impurity scattering 
theory [19]. In Figs. 8(a)-(c) we show our calculated 
p{T) for three representative carrier densities (all for the 
Si-22 sample of rcf. [22]) n = 5.7, 10.5, 21.3x10"cto-2 
including realistic surface roughness and interface (ox- 
ide) charged impurity scattering. At higher densities the 
quantitative agreement between the realistic calculations 
(solid lines) including both charged impurity and sur- 
face roughness scattering and the experimental data is 
very good whereas at the lower density (5.7x lO^^cm^^) 
the agreement is systematically poorer at higher tem- 
peratures where the theory consistently underestimates 
the experimental temperature dependence. However, if 
we use the actual density dependent effective mass mea- 
sured in the experiment rather than the band mass, we 
obtain much better agreement between our result and 
the experimental data at low densities. In Fig. 8(c) we 
show our calculated resistivity (thick solid line) using the 
measured effective mass taken from ref. [24]. We note, 
however, that the effective mass renormalization (and its 
effect on transport) is a rather subtle issue, and, in spite 
of the excellent quantitative agreement obtained in Fig. 



8(c) using the experimentally measured density depen- 
dent effective mass, it is unclear whether the band mass 
or the renormalized mass should be used in our RPA- 
Boltzmann transport theory. 

The quantitative agreement shown in Fig. 8 can be fur- 
ther improved by including hulk charged impurity scatter- 
ing due to random ionized impurity centers in the semi- 
conductor material itself, whose relative importance in- 
creases at lower densities as the 2D electrons reside on 
the average more deeply inside Si (and away from the 
interface). If we include in addition scattering by neutral 
impurities (with short-range bare impurity potential in 
contrast to the long-range bare impurity potential aris- 
ing from the charged impurities), we can essentially ob- 
tain quantitative agreement between theory and experi- 
ment by suitably adjusting the relative impurity densities 
among oxide charged impurities, bulk ionized impurities, 
and neutral short-range impurities (in addition, of course, 
to short-range siuface roughness scattering). But such a 
parametrized quantitative agreement between theory and 
experiment is essentially a device simulation exercise and 
is completely meaningless from the perspective of the fun- 
damental physics of 2D carrier systems, since it does not 
tell us anything more than what the results of Figs. 7 and 
8 already tell us. Th(> important point to realize here is 
that even the qualitatively valid results of Fig. 7 already 
show the same level of quantitative agreement between 
the RPA-Boltzmann theory and the Si-22 data of ref. 
[22] as the existing agreement among the experimental 
data from different Si MOSFET samples at same densi- 
ties and temperatures — this can be easily seen just by 
comparing to the data in Fig. 1(a) of ref. [22] correspond- 
ing to the Si-22 sample with those in Fig. 1(b), (c) of ref. 
[22] corresponding to the Si- 15 sample (in fact, sample 
to sample variations in the observed p{T, n) in various Si 
MOSFETs are typically much larger than the quantita- 
tive agreement we find without adjusting any parameters 
in Fig. 8). Our results shown in Figs. 7 and 8 also re- 
inforce the real danger of insisting just on "explaining" 
theoretically the low-temperature slope of the resistiv- 
ity (which automatically assumes a linear temperature 
dependence in both resistivity and conductivity) as has 
been fashionable in the recent experimental work [11,22] 
motivated by the asymptotic low temperature analysis of 
the interaction theory [21]. As Figs. 7 and 8 explicitly 
demonstrate, an agreement between theory and experi- 
ment on dp/dT in the T ^ Q limit is absolutely no guar- 
antee for a quantitative agreement between theory and 
experiment in a reasonable temperature range. Finally, 
we point out that the unrealistic ideal 2D approxima- 
tion and bare zero-range disorder potential misleadingly 
provides "better" quantitative agreement between theory 
and experiment, which of course does not mean anything 
about the physics of 2D carrier systems. We also note in 
this context that the quantitative agreement between the 
RPA-Boltzmann theory and the low density experimen- 
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tal data can be substantially improved [12] by assuming 
an effective lower density for the free carriers (participat- 
ing in "metallic" transport) than the nominal 2D carrier 
density, which is equivalent to invoking a lower effective 
Fermi temperature for the system (this could also arise 
from an enhanced effective mass at lower densities as has 
been experimentally observed), thus effectively enhanc- 
ing the theoretical temperature dependence. 

It is instructive to separate on the possible reasons 
for the systematic deviation of the experimental p{T) 
from the RPA-Boltzmann theory at lower densities (< 
5xl0^^cm~^) and moderate temperatures (> 2K). There 
arc two general possibilities: limitation or shortcomings 
of the RPA-Boltzmann theoretical scheme and the par- 
ticipation of scattering mechanisms left out of our model. 
The RPA-Boltzmann transport theory is obviously quite 
approximate since it is fundamentally semiclassical in na- 
ture leaving out all interaction, localization, and quan- 
tum interference effects. It is, however, difficult to un- 
derstand why these effects (i.e., localization, interaction, 
interference), which are left out of the RPA-Boltzmann 
transport theory, would become quantitatively more im- 
portant at higher temperatures. In fact, quantum effects 
should be progressively weaker as T/Tp increases, and 
therefore the systematic underestimation of the experi- 
mental p{T) in our RPA-Boltzmann theory is very un- 
likely to arise from the theoretical approximations of our 
model. We believe that the systematic apparent disagree- 
ment between experiment and theory at higher tempera- 
tures and lower densities in Fig. 8, if real, must arise from 
a missing scattering mechanism (neglected in our model) 
which takes on significance at higher temperatures. Such 
additional (high temperature) scattering may be due to 
intervalley scattering in Si and/or surface/oxide phonon 
modes not included in our theory. 

In discussing the theoretical approximations of our 
model it is important to emphasize, particularly in the 
context of the results shown in Figs. 7 and 8 for com- 
parison with the data of Ref . [22] , that the leading-order 
single-site scattering approximation used in our theory 
is quite accurate since the typical values of ^kpV (where 
kp and I are respectively the Fermi wave vector and the 
carrier mean free path) are rather large for the results 
depicted in Figs. 7 and 8. It is easy to show that for 
a 2D system kpl = ^IQv, where a = (j/{e^ /h) is the 
dimensionless conductance of the 2D system. For the 
results of Ref. [22], as shown in Figs. 7 and 8 of our 
paper, kpl is typically 100 or higher, making the single- 
site approximation well-valid. Thus, we do not believe 
that higher-order impurity scattering effects (or localiza- 
tion effects) are important in the regime of density and 
temperature covered in our Figs. 7 and 8. Obviously, at 
lower densities, near the critical density for the metal- 
insulator transition (where kpl < 1), higher-order impu- 
rity scattering effects become very important, and our 
simple RPA-Boltzmann theory breaks down. 



VI. DISCUSSION AND CONCLUSION 

We have carried out a critical comparison between the 
hiW p{T) calculated within the RPA-Boltzmann theory 
and its asymptotic (T/Tp — s- 0) linear-T approximation, 
finding that the linear approximation strictly applies only 
at very low temperatures, typically T/Tp < 0.05, at least 
within the RPA screening theory. In addition, we have 
provided a qualitative explanation, based on our approx- 
imate finding Ap{T,n) » Ap(f,(/o) where t = T/Tp 
and go = (iTp/'^kp, for the relative strength of metal- 
licity in various 2D systems (showing in the process that 
the dimensionless screening parameter and not the 
dimensionless interaction parameter rg, provides a bet- 
ter zeroth-order qualitative description for the metallicity 
trend in various materials). We have also carried out a 
detailed analysis of p{T,n) in the strictly 2D limit {i.e., 
with the subband form factor / = 1) and using short- 
range bare disorder. The strictly 2D results are in general 
quantitatively incorrect showing a more prominent lin- 
ear low-T regime and manifesting much stronger metal- 
licity than the realistic quasi-2D results. Similarly, the 
short-range bare disorder manifests stronger temperature 
dependence than the realistic charged impurity scatter- 
ing case although they both have the same leading-order 
temperature dependence. We have also carried out our 
theoretical transport calculations (not presented in this 
paper) including local field corrections to 2D screening 
by going beyond RPA (which is exact only in the high 
density or in the high temperature limit). Our results 
with local field corrections are qualitatively very simi- 
lar to the RPA screening results presented in this paper 
except that the temperature dependence of p{T) is in 
general somewhat weaker. We restrict ourselves to pre- 
senting only RPA screening results because there is no 
unanimity in the literature about the best possible local 
field corrections and also because the temperature de- 
pendence of local field effects are in general unknown. 
We point out that recent experiments attempting to ver- 
ify the interaction theory [11,22,23] have produced con- 
flicting results mainly because a clear linear-T regime in 
conductivity satisfying Tjj <^T <^Tp seems not to exist 
in low density "metallic" 2D systems (as wc explicitly 
show for the screening theory results in this paper) in 
any experimentally accessible range of the low-T data. 
One should therefore be extremely careful in applying 
any leading-order asymptotic temperature expansion to 
the 2D resistivity in analyzing experimental data. 

This caution in comparing 2D transport data with the 
interaction theory [21] is particularly warranted in view 
of the minimal necessary (but by no means sufficient) 
conditions that must be satisfied for such a comparison 
to be meaningful: (1) The 'weak' temperature constraint 
Td <C T <C Tp must be obeyed so as to be in the low tem- 
perature ballistic limit and phonon scattering must be 
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negligible, (2) a clear linear temperature dependence in 
the low temperature conductivity (noi resistivity) a{T) = 
ao[^+C-l{T/TF)], where the slope Ci depends on density, 
must be observed over a reasonable (a decade or more) 
range of temperature satisfying the ballistic constraint 
Tjj <^ T <^ Tp; (3) the actual temperature dependent 
conductivity correction, Aa/aQ = \(y(T) — (7o|/(To, must 
be very small (Aa/ao <C 1) for the leading-order interac- 
tion theory to be applicable. This set of necessary con- 
ditions is sufficiently restrictive so that no 2D transport 
experiment actually satisfies all three conditions except 
at very high carrier densities where the RPA-Boltzmann 
theory gives quantitatively accurate results [27]. 

Before concluding we discuss the experimental rele- 
vance of our (go, t) scaling prediction and the theoretical 
relevance (or validity) of our RPA screening approxima- 
tion. A cursory examination of the available 2D trans- 
port data shows that our predicted scaling of Ap/ po with 
the parameters qq = qxp/'^kp and t = T/Tp works on a 
qualitative level as a zeroth order description for GaAs 
electrons and holes, and also for electrons in Si MOS but 
only at low values of T/Tp (< 0.2). At higher tempera- 
tures phonon scattering becomes significant [13,16,23] in 
2D GaAs systems whereas in Si MOS structures, where 
phonon effects should be negligible, the temperature de- 
pendence of resistivity at lower metallic densities be- 
comes stronger than our prediction for reasons unclear 
to us. We note that this systematic under-estimation of 
experimental temperature dependence of p{T, n) at lower 
(higher) values of n (T) in the RPA scattering theory is 
unlikely to be arising from the interaction effects con- 
sidered in Ref. [21] since the disagreement arises in the 
non-asymptotic regime of T/Tp > 0.2 where the experi- 
mental p{T, n) is manifestly non-linear in T/Tp making 
the theory of Ref. [21] inapplicable. We can speculate 
several possible reasons for this unusually strong tem- 
perature dependence of Si MOSFETs: (1) Somehow the 
effective Fermi temperature (and the Fermi wave vec- 
tor) in the system could be smaller than the nominal 
Fermi temperature (or the Fermi wave vector) obtained 
on the basis of (100) Si inversion layer band mass and 
measured carrier density (leading to enhanced values of 
50 and/or t) - for example, the effective mass could be 
larger due to renormalization or the effective free car- 
rier density smaller (e.g. due to trapping by interface 
defects); (2) there could be additional scattering mech- 
anisms (e.g. intersubband scattering between different 
valleys in Si, which would be enhanced at higher tem- 
peratures) not included in our theory; (3) RPA could be 
failing systematically at higher values of T/Tp (this is 
an unlikely scenario since RPA should become a better 
approximation at higher T/Tp values, and going beyond 
RPA including local field corrections does not help in this 
respect). 

We note that we have ignored in this work any damping 
correction to the RPA screening function arising from the 



impurity scattering effect. In particular, for T <^ To the 
temperature-induced thermal suppression of 2kp screen- 
ing, which is crucial in producing the strong tempera- 
ture dependence of resistivity [19] at low temperature, 
becomes completely ineffective since collisional broaden- 
ing or damping effects induced by impurity scattering 
already suppresses screening at q = 2k p. Therefore, for 
T <C Id we expect the 2D resistivity p{T) to become es- 
sentially independent (or a very weak sublinear function) 
of temperature. This is, in fact, precisely the experimen- 
tal observation — the strong (and often approximately 
linear) temperature dependence of the low temperature 
resistivity almost always shows a saturation at very low 
temperatures for T <C To- This damping or broadening 
induced screening suppression for T <^ Tjj has a more 
subtle effect also. At low carrier densities Tp is low and 
Tp) is typically high since scattering effect is strong at 
low density; therefore at very low densities, a "metallic" 
2D system may not manifest strong "metallicity" because 
the damping induced low temperature saturation of p{T) 
will be more important as Tn approaches Tp at low den- 
sities. This effect is also consistent with experimental ob- 
servations. We have discussed elsewhere [12] the damp- 
ing correction to screening in some details in the con- 
text of the 2D MIT phenomenon. In this article we have 
restricted ourselves mostly to discussing the "ballistic" 
T ^ Td regime of 2D metallicity, and as such, we have 
decided to ignore the damping correction. The other rea- 
son for ignoring the damping correction is that the precise 
Dingle temperature value Tjj to be used in the theory is 
unknown, and therefore it introduces an unknown free 
parameter which we wish to avoid. Also, our purpose 
of comparison with the interaction theory of Zala et al. 
[21] is not well-served by having the damping correction 
since the interaction theory has been explicitly developed 
for the ballistic regime. We do mention, however, that 
introducing a Dingle temperature induced damping cor- 
rection to the RPA screening function will, in general, 
reduce the overall temperature dependence of our cal- 
culated resistivity with the low-temperature {T <C Tp)) 
resistivity showing an approximate saturation behavior 
[12]. 

We emphasize that the screening theory produces ex- 
cellent qualitative agreement with the existing experi- 
mental data and provides a good zeroth order (go, t) 
scaling description, which is all one should expect from 
the simple Drude-Boltzmann RPA-screening theory for a 
strong-coupling (i.e. low carrier density) problem. The 
question of how valid Boltzmann theory is for under- 
standing 2D "metallicity" is a difficult question to an- 
swer. A serious problem in this respect is the fact that 
this is the only quantitative theory that exists in the 
literature for studying 2D metallicity, and therefore its 
validity can only be judged a posteriori through an em- 
pirical comparison with the experimental data. (There 
is simply no competing alternative quantitative theory 
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in this problem, making the discussion of the vahdity 
of the RPA-Boltzmann theory somewhat meaningless - 
we emphasize in this context that the interaction theory 
of ref. [21] is not an alternative theory, it is an exten- 
sion of the RPA theory to incorporate higher-order quan- 
tum interaction corrections within necessarily a highly re- 
strictive model; the interaction theory applies only when 
the temperature correction to conductivity is very small 
and linear.) We believe that the RPA screening descrip- 
tion should be qualitatively well valid in this problem as 
long as the dominant disorder in the 2D semiconductor 
systems arises from long-ranged random charged impu- 
rity scattering (as is the case here at low carrier densi- 
ties). This is because RPA is a physically motivated ap- 
proximation, leading to the screening of the long-ranged 
Coulomb scattering potential to a short-ranged screened 
disorder. The approximation becomes exact only in the 
limit of high density and high temperature, but the 2D 
mctallicity manifests itself only for finite values of T/Tp 
(i.e. for T/Tp extremely small p does not show strong 
T-dependence!) and as such, RPA should be a reason- 
able description. We believe that the exactness of RPA 
in the rs limit has been overemphasized in the lit- 
erature - RPA remains a qualitatively accurate descrip- 
tion of metallic systems even for large values of (par- 
ticularly at finite T/Tp) as long as there is no quan- 
tum phase transition to a non-Fermi liquid phase. The 
2D metal-insulator-transition (2D MIT) being essentially 
a "high-temperature" phenomenon, RPA, in our opin- 
ion, is a reasonable description. The fact that our cal- 
culated transport results change little by including lo- 
cal field corrections (beyond RPA) further reinforces the 
qualitative validity of RPA to this problem. It should, 
however, be emphasized that like any other (nonpertur- 
bative) uncontrolled approximation (e.g. DMFT, LDA) 
the quantitative accuracy of RPA is difficult to ascertain 
from a purely thcxjietical viewpoint. One advantage of 
the RPA-Boltzmann minimal transport model adopted 
in our work is that transport calculations can be carried 
out for arbitrary temperatures and carrier densities with 
concrete comparisons with experimental data. We have 
carried out one such comparison with a recent experi- 
ment [22] in this paper (sec. V), and earlier comparisons 
exist in the literature [14-16,19,23,26,27]. The general 
conclusion one can draw from these comparisons is that 
the RPA-Boltzmann theory is a reasonably successful ze- 
roth order model for 2D transport properties, providing 
a qualitative (and intuitively appealing) explanation for 
the strong 2D metallicity at low carrier densities. As 
a quantitative theory, however, it is not very successful 
at lower densities (which is not unexpected), and ad hoc 
theoretical refinements (e.g. assuming a lower effective 
carrier density or larger effective mass) may be needed 
for obtaining quantitative agreement with low density 
2D transport data. The systematic quantitative devi- 
ation of the leading-order RPA-Boltzmann theory from 



the experimental data on low-density 2D "metallic" sys- 
tems could arise from the large number of effects left out 
of this zeroth-order theory which take on significance as 
the carrier density is lowered, such as higher-order (be- 
yond screening) interaction corrections, higher-order im- 
purity scattering effects, possible localization corrections, 
and various additional scattering processes left out of the 
theory (e.g. intervalley scattering, bulk impurity scatter- 
ing, etc.) What is surprising is not that the zeroth-order 
RPA-Boltzmann theory becomes systematically quanti- 
tatively unreliable at lower carrier densities, but the 
fact that this minimal leading-order theory provides such 
an excellent qualitative description of the observed 2D 
metallicity (for example, by explaining the strong vari- 
ation in Ap(T, n) / po in various systems as arising from 
the difference in the qq, t values) at all densities and a 
good quantitative description at higher densities. This 
is particularly significant in view of the early suggestions 
made in the literature in the context of the 2D metal- 
insulator transition physics that the strong 2D metallic 
behavior must be arising from some exotic non- Fermi liq- 
uid ground state of the system at low carrier densities. 
It is now manifestly clear that the apparent 2D metal- 
lic behavior arises from standard Fermi liquid corrections 
involved in the interplay of interaction and disorder with 
screened effective disorder arising from the temperature 
dependent screening of random charged impurities in the 
system being the main qualitative source for the strong 
temperature dependent 2D resistivity. This high temper- 
ature consequences of the RPA-Boltzmann transport the- 
ory for 2D carrier systems are discussed in our compan- 
ion publication [26] with the current manuscript focusing 
entirely on the low- temperature transport properties. 
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APPENDIX A 

The carrier resistivity p in our theory is given by 

m 



P = 



ne'^{T) ' 



(Al) 



where m is the carrier effective mass, and the energy 
averaged transport relaxation time (r) is given in the 
Boltzmann theory by 



/<'eb(-1) 



where t{E) is the energy dependent relaxation time, and 
f{E) is the carrier (Fermi) distribution function. At T = 
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0, ,f{E) = 9{Ef — E) where Ep is the Fermi energy, and 
then (t) = t(Ep) giving the famihar result u = = 
ne^T{Ep)/m. The Fermi energy for a 2D system is given 
by Ep = imti? / [gytri) = ksTp, where a spin degeneracy 
of 2 has been assumed, and is the valley degeneracy 
of the 2D system {g^ = 1 for GaAs; 5^ = 2 for Si (100) 
MOSFETs). 

We calculate the impurity ensemble averaged relax- 
ation time t{E) due to elastic disorder scattering in the 
Born approximation: 



1 



27r 



T{Ek) h 



/OD 
dzNl"\z)\u^"\k--k';z)\'' 

x{l-cos9i,i,')5{Ek-Ek'), (A3) 



where E{k) = h^k^ /2m is the 2D carrier energy for 2D 
wave vector k; ^^k' is the scattering angle between car- 
rier scattering wave vectors k and k'; the delta function 
5{Ek — Ek') assures energy conservation for elastic scat- 
tering due to charged impurities where the screened scat- 
tering potential is denoted by u^"^ (q; z) with q = k — k' 
as the 2D scattering wave vector and z is the quanti- 
zation or the confinement direction normal to the 2D 
layer. The quantity n\°'\z) in Eq. (A3) denotes the 
3D charged impurity density (with the z dependence re- 
flecting a possible impurity distribution) of the a-th kind 
with a representing the various possible types of impu- 
rities which may be present in 2D semiconductor struc- 
tures. For example, a could denote charged impurities 
(or interface roughness) located at the interface (e.g. Si- 
SiOa interface for MOSFETs, GaAs-AlGaAs interface for 
GaAs heterostructures and quantum wells) or impurities 
in the 2D layer itself or remote charged impurities inside 
the insulator (which could nevertheless scatter the 2D 
carriers by virtue of the long-range nature of Coulomb 
scattering). It is, for example, known that in low density 
Si inversion layer's the dominant scattering arises from 
the charged impurities located at the Si-Si02 interface. 
In general, however, the distribution of scattering centers 
in 2D semiconductor systems is not known, and has to be 
inferred from a careful comparison between experimental 
transport data and theoretical calculations assuming var- 
ious kinds of scatters. Fortunately for our purpose, this 
lack of precise knowledge of the charged impurity distri- 
bution is not a crucial issue since, on the qualitative level 
of our interest, all of them give rise to strong tempera- 
ture dependence of the resistivity provided the carrier 
density is low enough - the temperature dependence is 
the strongest when the impurities are distributed in the 
bulk of the 2D layer (bulk impurities) and goes down 
as the impurities are located further (remote impurities) 
from the layer inside the insulator. In our calculation 
presented in this paper, we typically assumes only one 
kind of scatters parameterized by a single impurity den- 
sity to keep the number of parameters a minimum - this 
impurity density essentially sets the overall scale of resis- 



tivity in our results. We emphasize that we can obtain 
good qualitative agreement with almost all the existing 
2D MIT experimental data by choosing three different 
kinds of charged impurities (i.e. interface, remote, and 
bulk) parameterized by a few reasonable parameters, but 
we do not see much point in this data fitting-type en- 
deavor. 

The screened impurity potential '«*^")(q; z) = u'^"\q; z) 
is given by: 



U^-\q-z)^Wr'V.!^l{q-.z) 



(A4) 



where V^|^p is the bare potential due to a charged impu- 
rity and e{q) is the carrier dielectric function screening 
the impurity potential. The bare potential is given by 



27rZ(«)e2 
Rq 



(A5) 



where Z^"^ is the impurity charge strength, R is the back- 
ground (static) lattice dielectric constant, and Fimp is 
a form factor determined by the location of the impu- 
rity and the subband wavefunction yj{z) defining the 2D 
confinement. (We do not show the explicit form of the 
form-factor Fimp here except to note that it reduces to 
^imp = e"^'^ for a strictly 2D layer, where |</'(2;)P = i5(z), 
with d being the separation of the impurity from the 2D 
layer - iov d = 0, i.e. when the impurity is in the layer, 
-Pimp = 1 in this pure 2D limit as one would expect.) 

The finite wave vector dielectric screening function is 
written in the RPA as 



e{q) = l-v{q)n{q,T), 



(A6) 



where v{q) = V2D{q)f{q) is the effective bare electron- 
electron (Coulomb) interaction in the system with 
V2D{q) ~ lue? j {kci) being the 2D Fourier transform of 
the usual 3D Coulomb potential, e^/(Kr), and /(g) be- 
ing the Coulomb form factor arising from the subband 
wavef unctions '>l}{z): 



/(g) = ^ / dz\ dz'\^{zt\^[z')\^ 



— OO -J — DC 



-q{z+z'\ 



(A7) 



The second term in Eq. (A7) arises from the image 
charge effect due to Ki ^ Kg, where Ki and the 
lattice dielectric constants of the insulator and the semi- 
conductor respectively (with R. = {Ki + Ks)/2). We note 
that in the strict 2D limit, when |V'(-2)P = f{<l) = 1- 
The 2D irreducible finite-temperature (and finite wave 
vector) polarizability function n(g;T) is given by the 
noninteracting polarizability (the irreducible "bubble") 
function within RPA: 



n((z,T) 



, n(g:r = o,//) 

cosh'^ ^ A*') ' 



(A8) 
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where (3 = {ksT)-'^. In Eq. (A8) Il{q:T = 0,Ef) is 
the zero-temperature noninteracting static polarizabihty 
given by: 



U{q;T = 0,EF) = NF 



2fcf 



0{q - 2kF) 



(A9) 

where Nf — gv'mj'K is the density of states at Fermi en- 
ergy, and kF = {2im/ g^Y^'^ is the 2D Fermi wave vector. 
The chemical potential ^, in Eq. (A8) at finite tempera- 
ture T is given by 



II = -\n[-l + e-K^{l3EF) 



(AlO) 



We note that the integration in Eq. (A8) is over the 
dummy variable /z' which is unrelated to the real Fermi 
energy Ef of the system. For going beyond RPA one 
would rewrite the noninteracting polarizabihty function 
n(g, T) to a model interacting polarizabihty function 
Ilint (?) T) which is written as 



ni„t(«,r) 



n(g,T) 



l-v{q)G{q,T)n{q,Ty 



(All) 



where 11 is the noninteracting polarizabihty function de- 
scribed above and G(g, T) is a suitable local field correc- 
tion which approximately incorporates correction eflPects 
neglected in RPA. 

It is now straightforward to see that under the con- 
ditions of strict 2D approximation, no remote impurity 
scattering, only one kind of impurity scattering (i.e. only 
one value of Zi and Ni characterizing the impurity scat- 
tering strength), and no local field corrections the calcu- 
lated resistivity p(T, n) of the system expressed as the di- 
mensionless quantity p/ po where po = p{T = 0) depends 
only on the variables go = (lTF/2kF (with the 2D Tomas- 
Fermi wave vector qtf = 2gyme^ / {Rh^)) and t = T/Tf, 
where 



t oc T{gym)n~^ . 



(Al2a) 
(A12b) 



For our fully realistic calculations (as well as for the ex- 
perimental system), however, this scaling relation, i.e., an 
exclusive {qo,t) dependence of resistivity, is violated due 
to the quasi-2D nature of the system (i.e. |V'(-z)P ^{^)): 
the presence of the insulator (i.e. Ks 7^ Ki); various types 
of impurity distributions in the 2D layer, the interface, 
and in the insulator; local field corrections, etc. We note 
that our calculation assumes one subband occupancy, 
i.e. only the ground 2D subband i^{z) is considered in 
our work. At higher temperatures (and lower densities) 
other (excited) subbands may get occupied by carriers 
in which case one would have to carry out a multisub- 
band transport calculation including intersubband scat- 
tering processes. Such a multisubband generalization of 



the Drude-Boltzmann formalism given above is straight- 
forward, but the actual calculation of p{T) becomes ex- 
tremely complicate in this situation, and has only been 
attempted recently in one special case by us [25] where 
we considered intersubband scattering between spin-split 
subbands in the valance band of p-GaAs 2D hole systems. 
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